Breakdown of the N=0 Quantum Hall State in graphene: two insulating regimes 
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We studied the unusual Quantum Hall Effect (QHE) near the charge neutrality point (CNP) 
in high-mobility graphene sample for magnetic fields up to 18 T. We observe breakdown of the 
delocalized QHE transport and strong increase in resistivities p^x, \pxy \ with decreasing Landau level 
filling for V <2, where we identify two insulating regimes. For 1 ^ li'l ^ 1/2 we find an exponential 



increase of 



,ji(H H^) ^j^jjjjj ^jjg range up to several resistance quanta Rk, while the Hall 



effect gradually disappears, consistent with the Hall insulator (HI) with local transport. Then, 
at ~ 1/2 a cusp in pxx{H) followed by an onset of even faster growth indicates transition to a 
collective insulator (CI) state. The likely candidate for this state is a pinned Wigner crystal. 

PACS numbers: 73.43.-f 73.63.-b 71.70.Di 
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Graphene is a honeycomb monolayer of C atoms, which 
forms the most common carbon allotrope - graphite, but 
has only recently been isolated in the two-dimensional 
(2D) form [Ij. A bipartite honeycomb sp^-bonded lattice 
gives graphene its unique electronic structure. Two dis- 
persion sheets associated with the electrons belonging to 
two different sublattices form the filled tt and the empty 
TT* bands, which meet at two distinct isolated points (val- 
leys K and K'), yielding point-like Fermi surfaces. Low- 
energy electronic states in each valley have a linear 2D 
conical dispersion, e(p) — vpp, and in addition to 2D 
momentum have a 2D pseudo-spin quantum number ac- 
counting for the two-sublattice structure. Such quasi- 
particles are formally described by the Dirac equation 
for chiral massless fermions and have peculiar transport 
properties [3, 

In a magnetic field H perpendicular to the graphene 
layer the spectrum of Dirac quasiparticles is quantized 
into Landau levels (LL) with energies En = ±huJcVN. 
Plus and minus signs correspond to electrons and holes, 
respectively, uJc = vpy/^eH/ (he) is the "cyclotron fre- 
quency" for Dirac fermions, and N = n' + 1/2 zL 1/2 
where n' — 0,1,2,... enumerates orbital wave functions 
and ±1/2 are pseudo-spin eigenvalues. For each En there 
are four states, corresponding to different spin and val- 
ley indices. In the presence of Zeeman interaction and 
inter-valley scattering, these states might split further, 
as illustrated in Fig. |4] (d). Unlike the case of Landau 
quantization for non-relativistic massive electrons, where 
En = huJc{N -f- 1/2) and lOc — eH/{mc), in graphene 
there is a field-independent level at i? = for = 0. At 
the CNP (in undoped graphene) this level is half-filled, 
being equally shared between particles and holes. Ex- 
perimentally, such peculiar LL structure is manifested in 
the unusual QHE observed in graphene, 0,01; where Hall 
conductivity is quantized as axy — 4(Z + 1/2) /Rk, I is an 
integer and Rk — h/e^ is the resistance quantum. The 



nature of electronic states on the A^ = level remains un- 
clear and has recently become the focus of considerable 
attention d, H, 0, i, H • 

Here we report an experimental investigation of charge 
transport in a high-mobility graphene sample for low car- 
rier density n near the CNP and in magnetic fields up to 
18 T, that is, in the QHE regime near the A^ = Lan- 
dau level, whereprevious studies have yielded conflicting 
results d, 0, [1, H In particular, some studies have 
found finite longitudinal resistance p^x Rk near the 
CNP even for high magnetic fields, where the plateau at 
Pxy = Rk/2 around v — u^q/H = ±2 (<I>o is the flux 
quantum) corresponding to either particle or hole filling 
of the two lowest A^ LL is well developed 0, 0, S] ■ 
Others, reported pxx ^ Rk, in the Mfi range, indicating 
an insulating A^ = state at high fields In Ref. [13 

the Pxx{H) divergence with H was analyzed as an ad-hoc 
Kosterlitz-Thouless transition and associated with a crit- 
ical field He (rather than a filling Vc), which was found 
to be sample-dependent. The discrepancies between dif- 
ferent measurements could be somewhat reconciled by 
the fact that cleaner samples show stronger divergence 
of Pxx{H) with increasing H in the A^ = state [lo| . 
Hence, sample quality appears crucial for understanding 
the physics of the A^ = LL state in graphene. 

We have studied a monolayer graphene sample pre- 
pared by mechanical exfoliation of ZYA grade highly 
oriented pyrolytic graphite (HOPG) on a Si/Si02 sub- 
strate using the standard procedure described in [l[. 
Transport properties summarized in Fig. [1] (a,b) reveal 
high quality of our sample. Magnetoresistance measure- 
ments were performed using the low frequency (17.777 
Hz) lock-in technique with a driving current / = 10 nA 
and the 18 T superconducting magnet at the National 
High Magnetic Field Laboratory (NHMFL). The follow- 
up measurements in fields up to 6 T were performed at 
Purdue University. The sample mobility between the 
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FIG. 1: Longitudinal resistivity pxx of our sample in zero 
magnetic field as a function of the gate voltage Vg. (a) Solid 
lines are p^x at T = 5 K for increasing (thick) and decreas- 
ing (thin) Vg^ as indicated by arrows. Small hysteresis results 
from mobile charges in Si/Si02 substrate. The inset shows 
the same data together with the results of similar follow- 
up measurements at T = 2.5 K and T=150 K. Dotted lines 
are fits to p^x = + (epn)^. (b) Drude mobility, 

p = l/{nepxx), Q. The dashed line shows p = 2.82(6) • 10* 
cm^/V/s resulting from the fit. Filled and open symbols cor- 
respond to different Vg sweep directions shown by arrows. 
Triangles are Hall mobilities obtained from fits to the data in 
Fig. [U^d). The inset is the optical image of our sample. 

two measurements remained unchanged within the error, 
/i = 2.82(6) • 10'' cmVV/s, while the CNP has shifted 
by sa 2 V. At T = 150 K, p decreased by less than 10%, 
to p{150K) w 2.6 • 10'' cm^/V/s. A slight electron-hole 
asymmetry of the resistance in Fig. [l]is typical of devices 
with invasive contacts (see inset in Fig. [Tl (b)) and is ex- 
plained by the work function difference between graphene 
and the Au/Cr electrodes in our device [l^ . 

Fig. [5] summarizes the QHE in our sample for three 
charge carrier densities near the CNP, ?i = 1-10^^, 3.2- 
10" and 4.7 • 10" cm"^. The latter were refined by 
fitting the low-field linear part of the Hall resistivity to 
Pxy{H) = HI [nee), Fig. ^ (a). This yields dn/dVg = 
—8.07(9) • 10^° cm~^/V. Clear plateaux corresponding to 
half the resistance quantum in p^y and zero pxx develop 
for all three carrier densities upon approaching the LL 
filhng \v\ = 2, Fig. [2] (a) - (c). This is a hallmark of the 
QHE in graphene, resulting from only two of the total 
four iV = LL being available to either electrons, or 
holes, Fig. 2] (d). Hence, two conducting channels and 
Pxy = ±Rk/'^ plateaux. The developed QHE regime is 
also manifested by the plateau at 7r/2 in the Hall angle, 
^'Haii = ata.n{pxyl Pxx), Fig. [2] (d). Fitting the low-field 
linear part of ^Haii to the result of the Boltzman transport 
theory, tan(0Haii) — (m/c)^^, allows for an alternative 
refinement of the sample mobility p (dash-dotted lines in 
Fig. [11(d)). Thus obtained Hall mobility agrees very well 
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FIG. 2: QHE in our sample for different gate offsets Vg — Vd 
from the Dirac point Vd, i. e. for different carrier densities n. 
(a) pxy, (b) Pxx, (c) Uxy, aud (d) HaU angle, O^aii- u = ±2 fill- 
ing for each n is shown by the corresponding arrow. Dash-dot 
lines in (a) are linear fits of the low-field part of pxy{H) used 
to refine the HaU constant and obtain n = H /{ecpxy{H)). 
Curves in (b) are shifted for clarity, the zero level for each is 
given by the broken horizontal line. Horizontal dotted lines 
in (a) and (c) show the resistance and the conductivity quan- 
tization in graphene. 

with the Drude field mobility. Fig. [T] (b). 

When LL filling decreases far enough past = 2 
with the increasing magnetic field, the QHE resistance 
quantization breaks down and the system enters a resis- 
tive state, where Pxx{H) > and fluctuates widely with 
H. This fluctuating behavior of Pxx{H) and pxy{H) in 
Fig. d] (a), (b) is reproducible. In fact, the H > parts 
of the curves for n = 4.7 • 10^^ cm~^ overlay results of 
three different held sweep measurements, which essen- 
tially coincide. While v ^ 1 for this carrier density 
in our field range and Pxx.xy{H) increase only moder- 
ately, much lower = LL fillings, v <^ 1, are achiev- 
able for n = 10^^ cm^^. Here we observe dramatic in- 
crease of Pxx,xy{H), which is detailed in Fig. [31 In terms 
of the Hall conductivity it appears as a "zero plateau" 
state with axy — around v ~ 0, Fig. [D (c). It is im- 
portant to notice, however, that the Onsager symmetry 
Pxy{H) — —pxy{—H) is violated in this state, indicating 
breakdown of the dissipation-less QHE transport via de- 
localized states, typical of a Hall insulator (HI) [l3'|. This 
behavior is also reflected by the Hall angle dependence 
in (d), which deviates from the |0Haii(^^)| = 7r/2 QHE 
plateau and tends towards zero in high magnetic field. 

Fig. [3] shows the full range of Pxx(H) and pxy{H) vari- 
ation for n = 10^^ cm~^, whose initial parts are also in- 
cluded in Fig. [21 Careful examination of the data allows 
identifying distinct transport regimes. First, for some u 
in the range 1 < < 2 the QHE resistance quantization 
breaks down and a resistive state forms. The Hall an- 
gle becomes |6'Haii| < 7i'/2, but the Onsager relation still 
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FIG. 3: Breakdown of QHE in graphene for n ~ 10 cm~ 
and T = 0.25 K. Top scale shows magnetic field H correspond- 
ing to the inverse LL filling 1/u on the bottom. Filling per val- 
ley is Uk,k' ~ (a) pxx{H) (dark) and pxy{H) (light); (b) 
dpxx{H) / dH . Dramatic increase in pxx{H) beyond a cusp at 
~ ARk and a step- like anomaly in dpxx{H) / dH indicate tran- 
sition to a collective insulator state at ~ ±1/2. (c) blow-up 
of the initial growth in pxx,xy{H) and (d) the sum (dark) 
and the difference (light) of pxy for two opposite directions 
of magnetic field show disappearance of the Hall component 
following breakdown of |!/| = 2 QHE state with increasing H 
(see also Fig. [2Ia,b)). Dash-dotted curves in (c), (d) nearly 
coincident with the data in the HI phase at 1 > \v\ > 1/2 
are fits to an exponential increase, pa{H) = ae''^'^'"^"' [l5l |. 
Dotted curve in (c) is the difference, pxx{H) — pa{H). 

holds, as shown by the sum of Hall resistivities for two 
opposite magnetic field directions, pxy{H)+ pxy{—H) « 
in Fig. [3] (d) . Such behavior could be a sign of the full 
splitting of = level shown in Fig. [4] (d) and of a de- 
veloping v ^ 1 plateau, or of the unusual resistive Hall 
metal (HM) state at TV = LL considered in d, [li| . 

In the 1/2 < < 1 regime, a field-symmetric compo- 
nent breaking the Onsager relation appears in transverse 
resistivity pxy The Hall component pxy{H) + pxy{~H) 
decreases, while Pxx{H) and pxy(H) steadily increase up 
to several resistance quanta. Both Pxx{H) and the sym- 
metric part, pxy{H)+pxy{—H), follow the same exponen- 
tial dependence, Pa{H) — ae'^'^l^l^^'^-', shown by dash- 
dot lines in Fig. 0] (c,d). Independent fits in panels (c) 
and (d) give consistent He = 4.6(2) T, corresponding to 

w 1, or a filling vk.K' ~ v 12 k, 0.5 per valley. 

Such behavior can be understood as a Hall insulator 
(HI) [gI, [31 resulting from the Zeeman splitting of the 

= LL, Fig. 3] (d). In this case the delocalized Quan- 
tum Hall state and the related iV = mobility edge shift 
to a finite energy Eq k, gpsH (g is the Lande factor 
and pb is the Bohr's magneton). For filling factors be- 
low z/ « 1 the chemical potential falls below Eq and the 
electrons are localized. Thermally activated Hall trans- 
port via delocalized states is possible for T > 0, but it 
vanishes with the increasing field-induced spin splitting 
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FIG. 4: Breakdown of the Quantum Hall state in graphene 
at Ji" = 18 T as a function of gate voltage (carrier concentra- 
tion) at T = 0.25 K. (a) Onset of a strong increase of pxx and 
\pxy \ beyond ~ ARk near filling p — ±1/2 (shown by arrows). 
Blow-up of Pxy in (b) and axy in (c) show a well-developed 
QHE plateau around v = —2 and a hint of plateaux develop- 
ing at = 1 at Rk and 1/Rk, respectively. The apparent 
V = Q plateau in (c) corresponds to a bulk insulator with zero 
Hall angle and no Hall effect, (d) Schematics of Landau lev- 
els in graphene, including Zeeman splitting (up/down arrows) 
and K, K' valley splitting (bars) . 



of the N = level. At T w the system is insulating, 
the transport occurs via hopping between localized states 
and is dominated by the mesoscopic conductance fluctu- 
ations in the sample. The transverse voltage drop results 
from the sample and the leads average asymmetry. The 
I — V curve is expected to be strongly non-linear, as it 
was indeed observed in Ref. [13, although authors there 
have interpreted this nonlinearity as resulting from sam- 
ple heating. A similar = insulating regime in a lower 
mobility sample was recently reported in Ref. 

Probably the most surprising and intriguing is the 
;/ < 1/2 regime, which is characterized by dramatic 
changes in behavior oi pxx{H) and pxy{H). Both Pxx{H) 
and the symmetrized pxy{H) + pxy{—H) deviate from the 
exponential growth describing the HI phase. The differ- 
ence Pxx{H) — Pa{H) shown by the dotted line in Fig. [3] 
(c) emphasizes a cusp in pxxiH), which is also visible in 
Fig. [3] (a), and is followed by an abrupt increase of re- 
sistance beyond pxx ~ ^Rk ■ The cusp-like singularity in 
Pxx{H) is also identified by the prominent jump in the 
derivative dpxx/dH at j/ « 1/2, Fig. [3] (b). The Hall 
component of the transverse resistivity given by the dif- 
ference pxy{H) — pxy{—H) also shows singular behavior, 
abruptly approaching zero at the same filling. Fig. [3] (d) , 
so that pxy{H) « pxy{—H) for i/ < 1/2. 

Simultaneous singular behavior of pxy{H) — pxy{—H) 
and dpxx/dH at ~ 1/2 indicate transition to a dif- 
ferent, more insulating state, which we identify as col- 
lective insulator (CI). A likely candidate for such a CI 
state is a pinned Wigner crystal (WC), where electrons 
are collectively, rather than individually localized. Our 
interpretation is based upon analogy with QHE in high- 
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mobility 2D electron gases (2DEG) in semiconductor het- 
erostructures, where an onset of strongly insulating be- 
havior at low filling factors, v < 1/4, has also been re- 
ported 3,12,3- This collective- insulating behavior in 
2DEG is understood as a 2D WC that can be pinned by 
arbitrarily small disorder at zero In our case 

of graphene transition at v k, 1/2 corresponds to a filling 
vk,k' ~ 1/4 per valley, roughly consistent with 2DEG 
results. 

Finally, we have also investigated the breakdown of 
the A*" = Quantum Hall state and the appearance of 
the insulating behavior in our graphene sample by vary- 
ing the gate voltage Vg in magnetic field = 18 T. The 
results shown in Fig. [?] are in agreement with those in 
Fig. [21 supporting the picture presented above. There 
is a well-developed QHE plateau around filling \v\ — 2. 
The quantization breaks down in the 1 ^ ^ 2 range 
and there is a hint of a plateau at \pxy\ — Rk devel- 
oping around ji'l « 1, Fig. [4] (b), (c). An onset of 
strong increase of pxx and \pxy\ beyond ~ Rk is seen 
at w 1/2. It can also be traced in the behavior of the 
derivative dpxy/dVg. The fact that such an abrupt onset 
occurs here and in Fig.[3]at very different \n\ and \B\ but 
similar — \nh/eH\ w 1/2 , suggests that it is driven 
by the correlation between electrons in graphene, again 
consistent with a transition to a bulk collective insulator 
at < 1/2. Due to the contact-induced electron-hole 
asymmetry, there is a slight shift of the curves in Fig. [4] 
with respect to the nominal n = 0, making such mea- 
surements by sweeping Vg less favorable compared to the 
sweep of magnetic field at a constant Vg. 

In summary, we have investigated the breakdown of the 
Quantum Hall Effect and the emergence of an insulating 
behavior in the N = Landau Level in a high-mobility 
single-layer graphene sample. The LL filling in the range 
\v\ < 2 is achieved by either increasing the magnetic 
field at a constant carrier density n, or by varying n{Vg) 
at = 18 T. Careful analysis of our data lead us to 
identify two different resistive regimes as a function of 
the decreasing LL filling i/ = n^o/H. First, the well- 
developed resistance quantization on the = 2 plateau 
breaks down and a dissipative state develops near the LL 
filling w 1. This can be understood as a HI resulting 
from the Zeeman splitting of the N = level, similar to 
that observed for higher N states 0, 0, El- For iV = 0, 
however, as a result there remain no delocalized states 
occupied by electrons (holes) at sufficiently low filling 
1^1 1 Fig- El (a), and the transport is local at T ^ 0. 
This observation agrees with recent findings reported in 
Ref. 15 and largely rules out the picture of a Hall metal 
with spin-polarized chiral currents Q, although not its 
refined version in [l3|, which effectively leads to a HI. 

Our most striking finding, though, is a well-defined 



onset of the dramatic resistance increase with decreasing 
filling at « 1/2. It is clearly defined by the anomalies in 
the field dependencies pxx{H) and pxy{H) and can be un- 
derstood as a transition to the collective bulk insulator 
state. Collective bulk insulating states have been pre- 
viously observed in 2DEG systems at low filling v and 
are comm only associated with pinned Wigner crystals 
[Hill, El, El US] - We suggest that a pinned WC is also 
likely candidate for the strongly insulating state which 
we have identified in our graphene sample for < 1/2. 
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